Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Characterization of the FKBP12-Encoding Genes in Aspergillus fumigatus.

Invasive aspergillosis, largely caused by Aspergillus fumigatus, is responsible for a growing number of deaths among immunosuppressed patients. Immunosuppressants such as FK506 (tacrolimus) that target calcineurin have shown promise for antifungal drug development. FK506-binding proteins (FKBPs) form a complex with calcineurin in the presence of FK506 (FKBP12-FK506) and inhibit calcineurin activity. Research on FKBPs in fungi is limited, and none of the FKBPs have been previously characterized in A. fumigatus. We identified four orthologous genes of FKBP12, the human FK506 binding partner, in A. fumigatus and designated them fkbp12-1, fkbp12-2, fkbp12-3, and fkbp12-4. Deletional analysis of the four genes revealed that the Δfkbp12-1 strain was resistant to FK506, indicating FKBP12-1 as the key mediator of FK506-binding to calcineurin. The endogenously expressed FKBP12-1-EGFP fusion protein localized to the cytoplasm and nuclei under normal growth conditions but also to the hyphal septa following FK506 treatment, revealing its interaction with calcineurin. The FKBP12-1-EGFP fusion protein didn't localize at the septa in the presence of FK506 in the cnaA deletion background, confirming its interaction with calcineurin. Testing of all deletion strains in the Galleria mellonella model of aspergillosis suggested that these proteins don't play an important role in virulence. While the Δfkbp12-2 and Δfkbp12-3 strains didn't show any discernable phenotype, the Δfkbp12-4 strain displayed slight growth defect under normal growth conditions and inhibition of the caspofungin-mediated "paradoxical growth effect" at higher concentrations of the antifungal caspofungin. Together, these results indicate that while only FKBP12-1 is the bona fide binding partner of FK506, leading to the inhibition of calcineurin in A. fumigatus, FKBP12-4 may play a role in basal growth and the caspofungin-mediated paradoxical growth response. Exploitation of differences between A. fumigatus FKBP12-1 and human FKBP12 will be critical for the generation of fungal-specific FK506 analogs to inhibit fungal calcineurin and treat invasive fungal disease.

Biochemical and Structural Studies on PrfA, the Transcriptional Regulator of Virulence in Listeria monocytogenes

Abstract

Listeria monocytogenes is a gram-positive soil saprophytic bacterium that is capable of causing fatal infection in humans. The main virulence regulator PrfA, a member of the Crp/FNR family of transcriptional regulators, activates the expression of essential proteins required for host cell invasion and cell-to-cell spread. The mechanism of PrfA activation and the identity of its small molecule coactivator have remained a mystery for more than 20 years, but it is hypothesized that PrfA shares mechanistic similarity to the E. coli cAMP binding protein, Crp. Crp activates gene expression by binding cAMP, increasing the DNA binding affinity of the protein and causing a significant DNA bend that facilitates RNA polymerase binding and downstream gene activation. Our data suggests PrfA activates virulence protein expression through a mechanism distinct from the canonical Crp activation mechanism that involves a combination of cysteine residue reduction and glutathione (GSH) binding.

Listeria lacking glutathione synthase (ΔgshF) is avirulent in mice; however virulence is rescued when the bacterium expresses the constitutively active PrfA mutant G145S. Interestingly, Listeria expressing a PrfA mutant in which its four cysteines are mutated to alanine (Quad PrfA), demonstrate a 30-fold decrease in virulence. The Quad and ΔgshF double mutant strains are avirulent. DNA-binding affinity, measured through fluorescence polarization assays, indicate reduction of the cysteine side chains is sufficient to allow PrfA to binds its physiological promoters Phly and PactA with low nanomolar affinity. Oxidized PrfA binds the promoters poorly.

Unexpectedly, Quad also binds promoter DNA with nanomolar affinity, suggesting that the cysteines play a role in transcription efficiency in addition to DNA binding. Both PrfA and Quad bind GSH at physiologically relevant and comparable affinities, however GSH did not affect DNA binding in either case. Thermal denaturation assays suggest that Quad and wild-type PrfA differ structurally upon binding GSH, which supports the in vivo difference in infection between the regulator and its mutant.

Structures of PrfA in complex with cognate DNA, determined through X-ray crystallography, further support the disparity between PrfA and Crp activation mechanisms as two structures of reduced PrfA bound to Phly (PrfA-Phly30 and PrfA-Phly24) suggest the DNA adopts a less bent DNA conformation when compared to Crp-cAMP- DNA. The structure of Quad-Phly30 confirms the DNA-binding data as the protein-DNA complex adopts the same overall conformation as PrfA-Phly.

From these results, we hypothesize a two-step activation mechanism wherein PrfA, oxidized upon cell entry and unable to bind DNA, is reduced upon its intracellular release and binds DNA, causing a slight bend in the promoter and small increase in transcription of PrfA-regulated genes. The structures of PrfA-Phly30 and PrfA-Phly24 likely visualize this intermediate complex. Increasing concentrations of GSH shift the protein to a (PrfA-GSH)-DNA complex which is fully active transcriptionally and is hypothesized to resemble closely the transcriptionally active structure of the cAMP-(Crp)-DNA complex. Thermal denaturation results suggest Quad PrfA is deficient in this second step, which explains the decrease in virulence and implicates the cysteine residues as critical for transcription efficiency. Further structural and biochemical studies are on-going to clarify this mechanism of activation.

Unfolded protein response genes regulated by CED-1 are required for Caenorhabditis elegans innate immunity.

The endoplasmic reticulum stress response, also known as the unfolded protein response (UPR), has been implicated in the normal physiology of immune defense and in several disorders, including diabetes, cancer, and neurodegenerative disease. Here, we show that the apoptotic receptor CED-1 and a network of PQN/ABU proteins involved in a noncanonical UPR response are required for proper defense to pathogen infection in Caenorhabditis elegans. A full-genome microarray analysis indicates that CED-1 functions to activate the expression of pqn/abu genes. We also show that ced-1 and pqn/abu genes are required for the survival of C. elegans exposed to live Salmonella enterica, and that overexpression of pqn/abu genes confers protection against pathogen-mediated killing. The results indicate that unfolded protein response genes, regulated in a CED-1-dependent manner, are involved in the C. elegans immune response to live bacteria.

Polymorphisms in the ACE and ADRB2 genes and risks of aging-associated phenotypes: the case of myocardial infarction.

Multiple functions of the beta2-adrenergic receptor (ADRB2) and angiotensin-converting enzyme (ACE) genes warrant studies of their associations with aging-related phenotypes. We focus on multimarker analyses and analyses of the effects of compound genotypes of two polymorphisms in the ADRB2 gene, rs1042713 and rs1042714, and 11 polymorphisms of the ACE gene, on the risk of such an aging-associated phenotype as myocardial infarction (MI). We used the data from a genotyped sample of the Framingham Heart Study Offspring (FHSO) cohort (n = 1500) followed for about 36 years with six examinations. The ADRB2 rs1042714 (C-->G) polymorphism and two moderately correlated (r(2) = 0.77) ACE polymorphisms, rs4363 (A-->G) and rs12449782 (A-->G), were significantly associated with risks of MI in this aging cohort in multimarker models. Predominantly linked ACE genotypes exhibited opposite effects on MI risks, e.g., the AA (rs12449782) genotype had a detrimental effect, whereas the predominantly linked AA (rs4363) genotype exhibited a protective effect. This trade-off occurs as a result of the opposite effects of rare compound genotypes of the ACE polymorphisms with a single dose of the AG heterozygote. This genetic trade-off is further augmented by the selective modulating effect of the rs1042714 ADRB2 polymorphism. The associations were not altered by adjustment for common MI risk factors. The results suggest that effects of single specific genetic variants of the ADRB2 and ACE genes on MI can be readily altered by gene-gene or/and gene-environmental interactions, especially in large heterogeneous samples. Multimarker genetic analyses should benefit studies of complex aging-associated phenotypes.

Recursive directional ligation by plasmid reconstruction allows rapid and seamless cloning of oligomeric genes.

This paper reports a new strategy, recursive directional ligation by plasmid reconstruction (PRe-RDL), to rapidly clone highly repetitive polypeptides of any sequence and specified length over a large range of molecular weights. In a single cycle of PRe-RDL, two halves of a parent plasmid, each containing a copy of an oligomer, are ligated together, thereby dimerizing the oligomer and reconstituting a functional plasmid. This process is carried out recursively to assemble an oligomeric gene with the desired number of repeats. PRe-RDL has several unique features that stem from the use of type IIs restriction endonucleases: first, PRe-RDL is a seamless cloning method that leaves no extraneous nucleotides at the ligation junction. Because it uses type IIs endonucleases to ligate the two halves of the plasmid, PRe-RDL also addresses the major limitation of RDL in that it abolishes any restriction on the gene sequence that can be oligomerized. The reconstitution of a functional plasmid only upon successful ligation in PRe-RDL also addresses two other limitations of RDL: the significant background from self-ligation of the vector observed in RDL, and the decreased efficiency of ligation due to nonproductive circularization of the insert. PRe-RDL can also be used to assemble genes that encode different sequences in a predetermined order to encode block copolymers or append leader and trailer peptide sequences to the oligomerized gene.

Ferrochelatase is a conserved downstream target of the blue light-sensing White collar complex in fungi.

Light is a universal signal perceived by organisms, including fungi, in which light regulates common and unique biological processes depending on the species. Previous research has established that conserved proteins, originally called White collar 1 and 2 from the ascomycete Neurospora crassa, regulate UV/blue light sensing. Homologous proteins function in distant relatives of N. crassa, including the basidiomycetes and zygomycetes, which diverged as long as a billion years ago. Here we conducted microarray experiments on the basidiomycete fungus Cryptococcus neoformans to identify light-regulated genes. Surprisingly, only a single gene was induced by light above the commonly used twofold threshold. This gene, HEM15, is predicted to encode a ferrochelatase that catalyses the final step in haem biosynthesis from highly photoreactive porphyrins. The C. neoformans gene complements a Saccharomyces cerevisiae hem15Delta strain and is essential for viability, and the Hem15 protein localizes to mitochondria, three lines of evidence that the gene encodes ferrochelatase. Regulation of HEM15 by light suggests a mechanism by which bwc1/bwc2 mutants are photosensitive and exhibit reduced virulence. We show that ferrochelatase is also light-regulated in a white collar-dependent fashion in N. crassa and the zygomycete Phycomyces blakesleeanus, indicating that ferrochelatase is an ancient target of photoregulation in the fungal kingdom.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Risk of ovarian cancer and inherited variants in relapse-associated genes.

BACKGROUND: We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk. METHODS AND FINDINGS: Women with and without invasive ovarian cancer (749 cases, 1,041 controls) were genotyped at 136 single nucleotide polymorphisms (SNPs) within 13 candidate genes. Risk was estimated for each SNP and for overall variation within each gene. At the gene-level, variation within MSL1 (male-specific lethal-1 homolog) was associated with risk of serous cancer (p = 0.03); haplotypes within PRPF31 (PRP31 pre-mRNA processing factor 31 homolog) were associated with risk of invasive disease (p = 0.03). MSL1 rs7211770 was associated with decreased risk of serous disease (OR 0.81, 95% CI 0.66-0.98; p = 0.03). SNPs in MFSD7, BTN3A3, ZNF200, PTPRS, and CCND1A were inversely associated with risk (p<0.05), and there was increased risk at HEXIM1 rs1053578 (p = 0.04, OR 1.40, 95% CI 1.02-1.91). CONCLUSIONS: Tumor studies can reveal novel genes worthy of follow-up for cancer susceptibility. Here, we found that inherited markers in the gene encoding MSL1, part of a complex that modifies the histone H4, may decrease risk of invasive serous ovarian cancer.

High-throughput identification of chemical inhibitors of E. coli Group 2 capsule biogenesis as anti-virulence agents.

Rising antibiotic resistance among Escherichia coli, the leading cause of urinary tract infections (UTIs), has placed a new focus on molecular pathogenesis studies, aiming to identify new therapeutic targets. Anti-virulence agents are attractive as chemotherapeutics to attenuate an organism during disease but not necessarily during benign commensalism, thus decreasing the stress on beneficial microbial communities and lessening the emergence of resistance. We and others have demonstrated that the K antigen capsule of E. coli is a preeminent virulence determinant during UTI and more invasive diseases. Components of assembly and export are highly conserved among the major K antigen capsular types associated with UTI-causing E. coli and are distinct from the capsule biogenesis machinery of many commensal E. coli, making these attractive therapeutic targets. We conducted a screen for anti-capsular small molecules and identified an agent designated "C7" that blocks the production of K1 and K5 capsules, unrelated polysaccharide types among the Group 2-3 capsules. Herein lies proof-of-concept that this screen may be implemented with larger chemical libraries to identify second-generation small-molecule inhibitors of capsule biogenesis. These inhibitors will lead to a better understanding of capsule biogenesis and may represent a new class of therapeutics.

Interaction of Cryptococcus neoformans Rim101 and protein kinase A regulates capsule.

Cryptococcus neoformans is a prevalent human fungal pathogen that must survive within various tissues in order to establish a human infection. We have identified the C. neoformans Rim101 transcription factor, a highly conserved pH-response regulator in many fungal species. The rim101 multiply sign in circle mutant strain displays growth defects similar to other fungal species in the presence of alkaline pH, increased salt concentrations, and iron limitation. However, the rim101 multiply sign in circle strain is also characterized by a striking defect in capsule, an important virulence-associated phenotype. This capsular defect is likely due to alterations in polysaccharide attachment to the cell surface, not in polysaccharide biosynthesis. In contrast to many other C. neoformans capsule-defective strains, the rim101 multiply sign in circle mutant is hypervirulent in animal models of cryptococcosis. Whereas Rim101 activation in other fungal species occurs through the conserved Rim pathway, we demonstrate that C. neoformans Rim101 is also activated by the cAMP/PKA pathway. We report here that C. neoformans uses PKA and the Rim pathway to regulate the localization, activation, and processing of the Rim101 transcription factor. We also demonstrate specific host-relevant activating conditions for Rim101 cleavage, showing that C. neoformans has co-opted conserved signaling pathways to respond to the specific niche within the infected host. These results establish a novel mechanism for Rim101 activation and the integration of two conserved signaling cascades in response to host environmental conditions.

Two genes on A/J chromosome 18 are associated with susceptibility to Staphylococcus aureus infection by combined microarray and QTL analyses.

Although it has recently been shown that A/J mice are highly susceptible to Staphylococcus aureus sepsis as compared to C57BL/6J, the specific genes responsible for this differential phenotype are unknown. Using chromosome substitution strains (CSS), we found that loci on chromosomes 8, 11, and 18 influence susceptibility to S. aureus sepsis in A/J mice. We then used two candidate gene selection strategies to identify genes on these three chromosomes associated with S. aureus susceptibility, and targeted genes identified by both gene selection strategies. First, we used whole genome transcription profiling to identify 191 (56 on chr. 8, 100 on chr. 11, and 35 on chr. 18) genes on our three chromosomes of interest that are differentially expressed between S. aureus-infected A/J and C57BL/6J. Second, we identified two significant quantitative trait loci (QTL) for survival post-infection on chr. 18 using N(2) backcross mice (F(1) [C18A]xC57BL/6J). Ten genes on chr. 18 (March3, Cep120, Chmp1b, Dcp2, Dtwd2, Isoc1, Lman1, Spire1, Tnfaip8, and Seh1l) mapped to the two significant QTL regions and were also identified by the expression array selection strategy. Using real-time PCR, 6 of these 10 genes (Chmp1b, Dtwd2, Isoc1, Lman1, Tnfaip8, and Seh1l) showed significantly different expression levels between S. aureus-infected A/J and C57BL/6J. For two (Tnfaip8 and Seh1l) of these 6 genes, siRNA-mediated knockdown of gene expression in S. aureus-challenged RAW264.7 macrophages induced significant changes in the cytokine response (IL-1 beta and GM-CSF) compared to negative controls. These cytokine response changes were consistent with those seen in S. aureus-challenged peritoneal macrophages from CSS 18 mice (which contain A/J chromosome 18 but are otherwise C57BL/6J), but not C57BL/6J mice. These findings suggest that two genes, Tnfaip8 and Seh1l, may contribute to susceptibility to S. aureus in A/J mice, and represent promising candidates for human genetic susceptibility studies.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.

Chromosomal organization of adrenergic receptor genes.

The adrenergic receptors (ARs) (subtypes alpha 1, alpha 2, beta 1, and beta 2) are a prototypic family of guanine nucleotide binding regulatory protein-coupled receptors that mediate the physiological effects of the hormone epinephrine and the neurotransmitter norepinephrine. We have previously assigned the genes for beta 2- and alpha 2-AR to human chromosomes 5 and 10, respectively. By Southern analysis of somatic cell hybrids and in situ chromosomal hybridization, we have now mapped the alpha 1-AR gene to chromosome 5q32----q34, the same position as beta 2-AR, and the beta 1-AR gene to chromosome 10q24----q26, the region where alpha 2-AR is located. In mouse, both alpha 2- and beta 1-AR genes were assigned to chromosome 19, and the alpha 1-AR locus was localized to chromosome 11. Pulsed field gel electrophoresis has shown that the alpha 1- and beta 2-AR genes in humans are within 300 kilobases (kb) and the distance between the alpha 2- and beta 1-AR genes is less than 225 kb. The proximity of these two pairs of AR genes and the sequence similarity that exists among all the ARs strongly suggest that they are evolutionarily related. Moreover, they likely arose from a common ancestral receptor gene and subsequently diverged through gene duplication and chromosomal duplication to perform their distinctive roles in mediating the physiological effects of catecholamines. The AR genes thus provide a paradigm for understanding the evolution of such structurally conserved yet functionally divergent families of receptor molecules.